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effect of the two now is considered. The phase spectrum 
IS a function defined as a + ^ - frequency, where a is the 
constant phase shift and b is the slope of the linear phase 
shift. Figure 1-16 shows the results of applying a 90- 
degree constant phase plus a linear phase component to 
the smusoids in Figure i-11. The zero-phase wavelet with 
the same amplitude spectrum as that in Figure MI was 
shifted in time by -0.2 s because of the linear phase shift, 
and converted to an antisymmetric form because of the 
constant 90-degree phase shift. 

Other variations in phase spectrum are shown in Figure 
1-17. The zero-phase wavelet (Figure M7a) can be 
modified to different shapes simply by changing the phase 
spectrum. It can be modified so much that it may not 
resemble the onginal wavelet shape as illustrated by the 
last example (Figure \-l7d). By keeping the amplitude 
spectrum unchanged, the wavelet shape can be changed 
by modifying the phase spectrum. 



1'2.3 Time-Domain Operations 



Consider a reflectivity sequence represented by the time 
senes (1, 0, 4). Also consider an impulsive source that 
causes an explosion at / = 0 with an amplitude of 1 The 
response of the reflectivity sequence to an impulse is 
called the impulse response. This physical process can be 
descnbed as follows: 



Time of 
Onset 



Reflectivity 
Sequence 



Source 



Response 



1 0 



1 



TTT 



One unit time later, suppose that the impulsive source 
generates an implosion with an amplitude of -i This 
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FIG. 1-17. The shape of a zero-phase wavelet (a) can be modified by introducing a nonzero-phase spectrum of any 
form as in (b), (c), (d). 



process is described as: 
Time of Reflectivity 

Onset Sequence Source Response 

~ TTl ^ 0 -i 

Note that the response is the reflectivity sequence scaled 
by impulse strength. Since a general source function is 
considered to be a sequence of explosive and implosive 
impulses, the individual impulse responses are added to 
obtain the combined response. This process is called 
linear superposition and is described in Table 1-1. 

In Table 1-1, the asterisk denotes convolution. The 
response of the reflectivity sequence (1, 0, ^) to the 
source wavelet (1, was obtained by convolving the 
two series. This is done computationally as shown in 
Table 1-2. 

A fixed array is set up from the reflectivity sequence . The 
source wavelet is reversed (folded) and moved (lagged) 
one sample at a time. At each lag, the elements that align 
are multiplied and the resulting products are summed. 
The mechanics of convolution are described in Table 1-3. 
The number of elements of output array cfA:j is given by m 
+ /? - 1, where m and n are the lengths of the operand 
array a(i) and the operator array b(j), respectively. 

When the roles of the arrays in Table 1-2 are inter- 
changed, the output array in Table 1-4 results. Note that 
the output response is identical to that in Table 1-2. 
Hence, it does not matter which array is fixed and which 
is moved as the convolution is performed. 

Seismic processing often requires measurement of the 
similarity or time alignment of two traces. Correlation is 
another time-domain operation that is used to make such 
measurements. Consider the following two wavelets: 



Wavelet 1: (2, 1, ~1, 0, 0) 
Wavelet 2: (0, 0,2, 1, -1) 

Although these wavelets are identical in shape, wavelet 2 
is shifted by two samples with respect to wavelet 1. The 
time lag at which they are most similar can be deter- 
mined. To do this, perform the operation on wavelet 1 as 
described in Table 1-3 without reversing wavelet 2 (omit 
Step 0). This is crosscorrelation and the result is shown in 
Table 1-5. Crosscorrelation measures how much two time 
series resemble each other. Crosscorrelation of a time 
series with itself is known as autocorrelation. 

From Table 1 -5 , note that maximum correlation occurs at 
lag -2. This suggests that if wavelet 2 were shifted two 
samples back in time, then these two wavelets would 
have maximum similarity. 

Table 1-6 shows the crosscorrelation values that result 
when the arrays are interchanged. This time the maxi- 
mum correlation occurs at lag 2. Thus, if wavelet 1 were 
shifted two samples forward in time, then these two 
wavelets would have maximum similarity. Also note that, 
unlike convolution, crosscorrelation is not commutative; 
i.e., the output depends on which array is fixed and 
which is moving (compare the results listed in Tables 1-5 
and 1-6). 

Table 1-7 shows the autocorrelation lags of wavelet 1. 
Note that maximum correlation occurs at zero lag, an 
important property of autocorrelation. Moreover, the 
autocorrelation function is symmetric. This is a property 
of real time series. Therefore, only one side of the 
autocorrelation needs to be computed. 

It is heuristically shown in Section 1 .2.4 that convolution 
in the time domain is equivalent to multiplication in the 
frequency domain (Bracewell, 1965). Since correlation is 
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Table 1-1. Linear superposition. 



Time of 
Onset 



Reflectivity 
Sequence 



Table 1-4. Convolution of reflectivity sequence (1, 0, |) with 
source wavelet (1, -|). 
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Expressed dififerently: 

(1,0, = (1 -ii-i) 
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Source 
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Response 
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Table 1-2. Convolution of source wavelet (1, -{) with the 
reflectivity sequence (1, 0, ^. 
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Reflectivity 
Sequence 



Output 
Response 



Table 1-5. Crosscorrelation of wavelet 1 with wavelet 2. 
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Table 1-3. Mechanics of the convolutlonal process. 

Fixed Array: 

fl(l) a{2) a{3) a(4) a{5) a(6) ad) a(8) 
Moving Array: 
b(\)b(2) bO) 

Given two arrays, a(i) and b(j): 
Step 0 = Reverse moving array bQ). 
Step 1 = Multiply in the vertical direction. 
Step 2 = Add the products and write as output 
point c(k). 

• Step 3 = Shift array b(j) one sample to the right 
and repeat Steps 1 and 2. 
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1 -I 

2 1 -1 

0 2 1-1 



1-10 0 



Output Lag 



0 2 1-1 
0 0 2 1-1 
0 0 2 1-1 
0 0 2 I -I 
0 0 2 1-1 
0 0 2 I -1 



-2 
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6 
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-2 
0 
0 
0 
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-2 
-1 
0 
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2 
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Table 1-6. Crosscorrelation of wavelet 2 with wavelet 1, 



2 1-1 0 0 
2 1-10 0 
2 1-10 
2 1 -I 
2 1 
2 



0 0 2 1-1 
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-1 0 

1 -1 



2 1 -1 
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0 
0 
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-2 
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6 
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2 1-10 0 -2 



-3 
-2 
-1 
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a convolution without reversing the moving array (see 
Table 1-3), a similar frequency-domain operation also 
applies to correlation. Figure M8 is a summary of 
frequency-domain descriptions of convolution and corre- 
lation. 

From Figure 1-18, note that both convolution and corre- 
lation produce an output with a spectral bandwidth that is 
common to both of the input series. The immediate 
example is the band-pass filtering process. Also note that 
phases are additive in case of convolution and subtractive 
in case of correlation (Bracewell, 1965). For autocorrela- 
tion, this implies that the output series is zero phase. This 
fact already was verified by the example in Table 1-7 
where it was shown that the autocorrelation is symmetric 
with respect to zero lag. 

As a measure of similarity, crosscorrelation is used 
widely at various stages of data processing. For instance, 
traces in a common-midpoint (CMP) gather are crosscor- 
related with a pilot trace to compute residual statics shifts 
(see Section 3.4). Again, the fundamental basis for com- 
puting velocity spectra is crosscorrelation. The building 
blocks of the Wiener filter (Section 2.6) are crosscorrela- 
tion of the desired output waveform with the input 
wavelet and autocorrelation of the input wavelet. 



Table 1-7. Autocorrelation of wavelet 1. 

2 1-1 0 0 Output Lag 

2 1 -1 0 0 0 -4 

2 1 -1 0 0 0 -3 

2 1-10 0 -2-2 

2 1-10 0 1-1 

2 1 -1 0 0 6 0 

2 1-10 0 II 

2 1-10 0 -2 2 

2 1-10 0 0 3 

2 1 -10 0 0 4 



One other important process is the vibroseis correlation. 
This involves crosscorrelation of a frequency-modulated 
source (sweep) signal with the recorded vibroseis trace. 
The sweep is a frequency-modulated vibroseis source 
input to the ground. The convolutional model for vibro- 
seis data is described in Section 2.7.7. Figure 1-19 shows 
a vibroseis sweep signal, a recorded common-source 
gather, and the correlated gather. The sweep length is 10 
s with a frequency band of 6 to 60 Hz. The 15-s 
uncorrelated vibroseis record yields a 5-s correlated 
record. Note that the eariy part of the uncorrelated 
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record contains low-frequency energy with increasingly 
higher frequencies at late times. This is because an 
upsweep (frequency increasing with time) signal was used 
in this data example. 

1.2.4 Frequency Filtering 

What happens to a wavelet when its amplitude spectrum 
is changed while its zero-phase character is preserved? 
To begin, consider the wavelet (summed trace 1) result- 
ing from superposition of two very low-frequency com- 
ponents in Figure 1-20. Then add increasingly higher 
frequency components to the Fourier synthesis (summed 
traces 2 through 5). Note that the wavelet in the time 
domain is compressed as the frequency bandwidth (the 
range of frequencies summed) is increased. Ultimately, if 
all the frequencies in the inverse Fourier transformation 
are included, then the resulting wavelet becomes a spike, 
as seen in Figure 1-21 (summed trace 6). Therefore, a 
spike is characterized as the in-phase synthesis of all 
frequencies from zero to the Nyquist. For all frequencies, 
the amplitude spectrum of a spike is unity, while its phase 
spectrum is zero. 

Figure 1 -22 shows five zero-phase wavelets, synthesized 
as shown in Figure 1-20. Note that all of them have band- 
limited amplitude spectra. A zero-phase band-limited 
wavelet can be used to filter a seismic trace. The output 
trace contains only those frequencies that make up the 
wavelet used in filtering. The time-domain representation 
of the wavelet is iht filter operator. The individual time 
samples of this operator are the filter coefficients. The 
process described here is zero-phase frequency filtering, 
since it does not modify the phase spectrum of the input 
trace, but merely band limits its amplitude spectrum. 

Frequency-domain filtering involves multiplying the am- 
plitude spectrum of the input seismic trace by that of the 
filter operator. The procedure is described in Figure 1-23. 
On the other hand, the filtering process in the time 
domain involves convolving the filter operator with the 
input time series. Figure 1-24 is a summary of the filter 
design and its time-domain application. The frequency- 
and time-domain formulations of the filtering process 
(Figures 1-23 and 1-24) are based on the following impor- 
tant concept in time series analysis: Convolution in the 
time domain is equivalent to multiplication in the fre- 
quency domain. Similarly, convolution in the frequency 
domain is equivalent to multiplication in the time domain 
(Bracewell, 1965). 

Frequency filtering can be in the form of band-pass, 
band-reject, high-pass (low-cut), or low-pass (high-cut) 
filters. All of these filters are based on the same principle: 
construct a zero-phase wavelet with an amplitude spec- 
trum that meets one of the four specifications. Band-pass 
filtering is used most because a seismic trace typically 
contains some low-frequency noise, such as ground roll, 
and some high-frequency ambient noise. The usable 
seismic reflection energy usually is confined to a band- 
width of approximately 10 to 70 Hz, with a dominant 
frequency around 30 Hz. 

Band-pass filtering is performed at various stages in data 
processing. If necessary, it can be performed before 
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FIG. 1-18. The frequency-domain description of convolution 
and correlation. 

deconvolution to suppress remaining ground roll energy 
and high-frequency ambient noise that otherwise would 
contaminate signal autocorrelation. Narrow band-pass 
filtering may be necessary before crosscorrelating traces 
in a CMP gather with a pilot trace for use in estimating 
residual statics shifts (Section 3.4). Band-pass filtering 
also can be performed before computing crosscorreia- 
tions during construction of the velocity spectrum for 
improved velocity picking (Section 3.3). Finally, it is a 
standard practice to apply a time-variant band-pass filter 
to stacked data (this section). 

Practical Aspects of Filter Design 

Application of a filter in the frequency or time domain 
(Figures 1-23 and 1-24) yields basically identical results. 
In practice, the time-domain approach is favored, since 
convolution involving a short array, such as a filter 
operator, is more economical than doing Fourier trans- 
forms. 

From Figure 1-22, the fundamental property of fre- 
quency filters can be stated as: The broader the band- 
width, the more compressed the filter operator: thus, 
fewer filter coefficients are required. This property also 
follows from the fundamental concept that the eff'ective 
time span of a time series is inversely proportional to its 
effective spectral bandwidth (Bracewell, 1965). 

In designing a band-pass filter, the goal is to pass a 
certain bandwidth with little or no modification, and to 
largely suppress the remaining part of the spectrum as 
much as practical. Initially, it appears that this goal can 
be met by defining the desired amplitude spectrum for the 
filter operator as follows: 

A f) = V' <f<fi 
[O, elsewhere, 

where/, and/2 are the cutoff frequencies. This is known 
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